Failure regimes in (1+) dimensions in fibrous materials 
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In this paper, we introduce a model for fracture in fibrous materials that takes into account the 
rupture height of the fibers, in contrast with previous models. Thus, we obtain the profile of the 
fracture and calculate its roughness, defined as the variance around the mean height. We investigate 
the relationship between the fracture roughness and the fracture toughness. 
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I. INTRODUCTION 

The fracture of a material under different loading con- 
ditions is a complicated phenomenon that continues to 
attract the attention of several physicists and engineers. 
The fracture follows a sequential process of nucleation, 
growth, and coalescence of numerous cracks. These 
cracks tend to create a rough profile which often can be 
described by self-affine scaling In the last decade, 
great theoretical and experimental efforts have been done 
trying to understand the process of formation and pro- 
pagation of cracks in materials. A great number of stud- 
ies have been devoted to this process from the point of 
view of statistical mechanics, that utilize concepts such 
as percolation, fractals, scaling law, etc. QJ^. A vari- 
ety of computational models for material fracture have 
been developed in recent years and provided interesting 
results However, the high degree of correlations 

between the constituents leads to a high computational 
cost. An alternative is to employ fiber bundle models in- 
troduced more than forty years ago ^ , in which bundles 
of unidirectional fibers form a system with low degree of 
correlations allowing the fracture process to be simulated 
in a large scale. 

The fracture of the material can occur through a big 
crack which percolates the sample. In this case we say 
that the fracture is catastrophic. When the fusion of 
small cracks cause the fracture of the material, it is 
called shredding |^ . In the catastrophic regime the frac- 
ture profile is reasonably smooth, while in the shredding 
regime it is very rough. A parameter of easy physical in- 
terpretation used to characterize fracture surfaces is the 
roughness. It is defined as the variance around of the 
mean height of the fracture profile, and measures the 
complexity of the crack path. Thus, the rougher is the 
fracture profile the harder is the crack path. The rough- 
ness has a direct relation to the fractal dimension, which 
characterizes the fractal character of the surface fracture. 
Therefore, a high fractal dimension indicates a very rough 
fracture profile [^]j7|. 

Another important parameter in the study of fracture, 
is the fracture toughness. This quantity measures the 



amount of energy that a material can absorb before it 
fractures. The fracture toughness is intimately related 
to the amount of cracks that appear in the material . 
Therefore, the larger the number of cracks in the sample 
more energy will be absorbed before catastrophic frac- 
ture occurs. There are several works in the literature 
which investigate the relationship between the fracture 
roughness and the fracture toughness 

In the present Rapid Communication, we simulated a 
model for fracture in fibrous materials which allows us to 
obtain the fracture profile of samples in contrast with pre- 
vious models |H,|l^-|l6|. We calculate the profile rough- 
ness and used the obtained values to define the transition 
between two fracture regimes: catastrophic and shred- 
ding. We also investigate the relationship between the 
fracture roughness and toughness for fibrous materials. 



II. MODEL 

Our model consists of a (l-l-l)d bundle of iVo parallel 
fibers all with the same elastic constant k. In order to 
simulate the height of the sample, the fibers are divided 
in 77 segments with the same length. The fiber bundle is 
fixed at both ends to two parallel plates. One plate is 
fixed, and in the other, a constant force F is applied, for 
example, by hanging a weight on it. This force is shared 
in the same amount, a, on each fiber of the bundle, which 
undergo the same linear deformation z = F/Nk, where 
N is the number of unbroken fibers. When the deforma- 
tion z reaches a critical value Zc, the failure probability 
of an isolated fiber is equal to 1. The failure probability 
of a fiber i is given by 



P,{S,t) = 



1 



■ exp 
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where rii is the number of unbroken neighboring fibers, 
5 = z/zc = F/Nkzc is the strain of the material, 
t ~ KbT/Ec is the normalized temperature, Kb is the 
Boltzmann constant, T is the absolute temperature, and 
Ec is the critical elastic energy. In this model, besides 
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finding the failure probability of a fiber, we have to in- 
dicate in which segment it breaks. This segment is ran- 
domly selected and the probability of the fiber to break 
in it is given by 



I) 



C 



(2) 



where rrij is a vector which indicates how many times a 
segment j broke and C — Tljij^j + !)• The form of the 
Eq. (|^) simulates a concentration of tension near to the 
region where the fiber bundle is more weak. 

At the beginning of the simulation, the bundle is sub- 
mitted to a initial strain given by 



F 



Nokzc 



(3) 



At each time step of the simulation we randomly choose 
a fiber of a set of Nq = qNo unbroken fibers. The number 
q represents a percentage of fibers and allows us to work 
with any system size. Then, using Eq. (|l|) we evaluate the 
fiber failure probability Pi and compare it with a random 
number r in the interval [0,1). If r < the fiber breaks. 
We then choose a segment j in the fiber and evaluate its 
probability (f) to break, using Eq. (||). If the probability 
(j) is higher than the random number p the fiber breaks 
in the chosen segment. If not, we analyze the neighbor- 
ing segments (i + 1) and (j — 1) and again, evaluate the 
probability cj). If the condition p < 4> does not hold to 
neither of the neighboring segments, we return to the ini- 
tial segment and test the condition p < (j) for a new value 
of p. This process continues until the condition p < (p is 
true. Once defined the segment where the fiber breaks, 
we begin to test all neighboring unbroken fibers. The 
first segment tested in the neighboring unbroken fibers 
is the segment where the previous fiber broke. The fail- 
ure probability Pi of these neighboring fibers increases 
due to the decreasing of rii and a cascade of breaking 
fibers may begin. This procedure describes the propa- 
gation of a crack through the fiber bundle, which occurs 
in the perpendicular direction to the applied force. The 
process of propagation stops when the test of the prob- 
ability does not allow rupture of any other fiber on the 
border of the crack or when the crack meets another al- 
ready formed crack. The same cascade propagation is 
attempted by choosing another fiber of the set Nq. After 
all the Nq fibers have been tested, the strain 5 is increased 
if some fibers have been broken. Since the force is fixed, 
the greater the number of broken fibers, the larger is the 
strain on the intact fibers and the higher is their failure 
probability. Then, another set of Nq unbroken fibers is 
chosen and the entire rupture process is restarted. The 
simulation terminates when all the fibers of the bundle 
are broken, i.e., when the bundle breaks apart. At this 
stage we consider the profile of the fracture and analyze 
its statistical properties. 



III. RESULTS 

Using the model described in the preceding section 
we simulated the fracture of a fibrous materials under 
a static force F. In these simulations we considered the 
elastic constant fc = 1, the critical deformation Zc = 1, 
the number of segments 77 = 100 and the normalized tem- 
perature t = 0.5. 

Figure |l] shows the fracture profile for a set of iVg = 600 
fibers and three different forces. Notice that, the lower 
the applied force the larger is the roughness of the frac- 
ture profile, i.e., more irregular is the crack path. For 
F — 4000 only one crack propagates in the material and 
for F — 2800 and F — 2000 more than one crack crosses 
the material. It can be seen in Fig. 0, that the rupture of 
the sample begins in different segments. This occurs be- 
cause we did not consider a deterministic starting notch 
in our simulations. 

Figure ||(a) shows the results obtained for the rough- 
ness W as a. function of the applied force F for different 
values of TVq. To measure the roughness of the fracture 
profile we used the method of the best linear least-square 
fitting described in Ref . [0 . In this plot we can observe 
that the roughness W decreases as the force F increases, 
and after a critical force F^ it stays constant. For each 
system size we have a characteristic value for Fc, which 
increases with the increase of the number of fibers iVg. 
This indicates that the greater the number of fibers in 
the bundle, the tougher is the sample. For F > Fc the 
fracture is catastrophic, i.e., the breakage of a fiber in- 
duces the rupture of the whole bundle. In this case the 
bundle breaks with only one crack. In this region the 
crack propagates in the material with high speed, leading 
to a quite rapid rupture process. For F < Fc the rupture 
of the bundle occurs due to the formation of small cracks, 
which weaken the bundle. Here, the crack speed tends to 
zero and a slow process of successive ruptures appears in 
the material. 

In Fig. ||(b), the x axis was converted to the initial 
strain Sq = F/Nq. The data are now found to collapse all 
on the same roughness-strain curve. From this curve we 
can find the critical initial strain Sqc, which does not de- 
pend on the system size. Our results indicate Sqc = 1-134. 
We may assume that, for an initial strain greater than 
6oc just one crack provokes the rupture of the material 
and that, for a initial strain below Sqc the rupture occurs 
through the fusion of small cracks. 

The critical value 6qc can be obtained analytically. The 
failure probability, Eq. can be written as 



P^iS,t) 



where r(t, 6) is defined as 



r(t, 6) — Sexp 



(4) 
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Using Eq. (|) we may observe that for T(t, 6) — 2.0 the 
breakage of any fiber induces the rupture of the bundle 
with just one crack. Then, we can assume that there is a 
critical value for r(i, S) = 2.0 that defines the transition 
between two regimes. In the first regime, a catastrophic 
fracture occurs in the first attempt to break the bundle, 
while in the second one the rupture of the bundle occurs 
due to the formation of small cracks, which weaken the 
bundle. From Eq. (||) we can show that the normalized 
temperature can be given by 

ln{Tc) - ln{5oc) 

So, for t = 0.5 and Tc{t,6) = 2.0, Eq. (|) will be valid 
only if Sqc « 1.134. 

Now we proceed to the evaluation of the fracture 
toughness kc- It can be defined by the work done to 
break the fiber bundle and is given by 

i^c^^n, (7) 

i 

where is the work done to break each fiber of the bun- 
dle. The work Ti is obtained by the following expression 

n = ^kzl (8) 

where Zi is the deformation of the fiber i. 

Figure || shows the log-log plot of the fracture tough- 
ness Kc versus the force F for different system sizes. This 
figure shows how the fracture toughness kc decreases with 
the increase of the force F, until it reaches a minimum 
value. This value is attained when a critical force Fc 
is applied to the system. Above of Fc the fiber bundle 
breaks catastrophically. It is known experimentally that 
a catastrophic break consumes little energy. As the force 
F decreases below Fc, the material can absorb more en- 
ergy before it fractures and the number of cracks in the 
material increases. Also, notice that the fracture tough- 
ness increases with the system size. Thus, the higher the 
number of fibers the more energy will be absorbed by 
fracture. 

Some works in the literature conjecture that there is a 
relationship between the fractal dimension and the frac- 
ture toughness |[^-p^, while others conjecture that such 
relationship does not exist |^,^. In this paper we in- 
vestigated the connection between the roughness W, di- 
rectly related to the fractal dimension, and the toughness 
fracture Kc- 

Figure ^ shows how the roughness W changes with 
the toughness Kc for several number of fibers Nq. In 
this figure we can see that the roughness W increases 
with the increase of the fracture toughness Kc, until it 
reaches a stationary state. We can interpret this figure 
in the following manner: with the increase of Kc, small 
cracks appear at all parts of the sample, making the frac- 
ture profile rougher. In the saturation region, we have a 



larger number of small cracks and they are totally un- 
correlated. The phenomenon of saturation constitutes a 
finite size effect and is related to the number of segments 
in which the fibers are divided. 



IV. CONCLUSION 

In conclusion, we studied a model for fracture in fi- 
brous materials in (1-1-1) dimensions which take into ac- 
count the rupture height of the fibers, in contrast with 
previous models. We obtained the fracture profile and 
evaluate its roughness and toughness. In this work we 
show that in the catastrophic regime the roughness of 
the fracture profile is reasonably smooth. In the shred- 
ding regime, in which slow cracks are formed in the ma- 
terial, the fracture profile is very rough. In this regime 
the energy necessary to break the material is higher than 
in the catastrophic regime. Our results indicate that the 
roughness W is related to the fracture toughness. We 
believe that the search for possible relationships between 
the roughness and the fracture toughness could stimu- 
late the further studies in order to check whether these 
relations are valid or not. 
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FIG. 1. Fracture profiles for three different forces. From 
left to right we have: F = 4000, F = 2800 and F = 2000 
arbitrary unity. In this simulation we used a total of 600 
fibers. 



1000 
800 
600 
400 

200 





1000 









»N0 = 2000 








am = 2200 




m □ \ % 




M.m = 2400 
vN0 = 2600 


- (a) 


* □ » & 

• □ J S 




1 >N0 = 2S00 




. □ J 

□ * 


t 
► 


\ *m = 3000 



FIG. 2. (a) Roughness W{F) as function of the applied 
force F for different number of fibers Nq. The data were av- 
eraged over 1000 samples, (b) The x-axis was converted to 
the initial strain 5o = F/Nq. 




1000 



1585 



logJF) 



2512 



3981 



FIG. 3. Fracture toughness Kc vs the applied force F for 

different number of fibers A^o. The data were averaged over 
1000 statistically independent samples. 
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FIG. 4. Roiiglmess W vs the fracture toughness Kc for dif- 
ferent number of fibers A'o. The data were averaged over 1000 
statistically independent samples. 
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